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ABSTRACT 

Three peaks and two dips have been detected in the power spectrum of the cosmic microwave 
background from the BOOMERANG experiment, at I ~ 210, 540, 840 and t ~ 420, 750, respec- 
tively. Using model- independent analyses, we find that all five features are statistically significant 
and we measure their location and amplitude. These are consistent with the adiabatic inflation- 
ary model. We also calculate the mean and variance of the peak and dip locations and amplitudes 
in a large 7-dimensional parameter space of such models, which gives good agreement with the 
model-independent estimates, and forecast where the next few peaks and dips should be found 
if the basic paradigm is correct. We test the robustness of our results by comparing Bayesian 
marginalization techniques on this space with likelihood maximization techniques applied to a sec- 
ond 7-dimensional cosmological parameter space, using an independent computational pipeline, 
and find excellent agreement: fi tot = 1.02±g;g| vs. 1.04±0.05, VL b h 2 = 0.022+°;™ 4 vs. 0.019+°;°°4, 
and n s — 0.96io'o9 vs - 0.90±0.08. The deviation in primordial spectral index n s is a consequence 
of the strong correlation with the optical depth. 

Subject headings: Cosmic Microwave Background Anisotropy, Cosmology 
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1. Introduction 

BOOMERANG has recently produced an im- 
proved power spectrum Ci of the cosmic microwave 
background (CMB) temperature anisotropy rang- 
ing from £ ~ 100 to I ~ 1000 (Netterfield et al. 
2001, hereafter B01). Three peaks are evident 
in the data, the first, at £ ~ 210 confirms the re- 
sults of previous analysis of a small subset of the 
BOOMERANG data set (de Bernardis et al 2000, 
hereafter BOO) as well as results from other exper- 
iments (see e.g. Miller et al. 1999; Mauskopf et al. 
2000; Hanany et al. 2000). Analysis of the bulk 
of the remaining data has improved the precision 
of the power spectrum and extended the coverage 
to £ ~ 1000, revealing for the first time a second 
peak at I ~ 540, and a third at I ~ 840. 

The results of two other experiments, released 
simultaneously with B01, are in good agreement 
with the BOOMERANG data. In addition to the 
first peak, the results from DASI (Halverson et al. 
2001) show a peak coincident with that seen in 
our data near £ ~ 540, and a rise in the spectrum 
toward high £ that is consistent with the leading 
edge of the peak seen in our data near £ <~ 840. 
The results from MAXIMA (Lee et al. 2001) are of 
lower precision, but are consistent with both the 
DASI and the BOOMERANG results. 

These detections are the first unambiguous con- 
firmation of the presence of acoustic oscillations 
in the primeval plasma before recombination (Pee- 
bles et al. 1970; Sunyaev & Zeldovich 1970), as ex- 
pected in the standard inflationary scenario (Bond 
& Efstathiou 1987). 

If the adiabatic cold dark matter (CDM) model 
with power-law initial perturbations describes our 
cosmogony, then the angular power spectrum of 
the CMB temperature anisotropy is a powerful 
tool to constrain cosmological parameters (see e.g. 
Kamionkowski & Kosowsky 1999, and references 
therein). 

In B01 a rigorous parameter extraction has 
been carried out with the methods of Lange et al. 
(2001), significantly improving the constraints ob- 
tained from previous CMB analyses (see e.g. Do- 
delson & Knox 2000; Melchiorri et al. 1999; Lange 
et al. 2001; Balbi et al. 2000; Tegmark & Zaldar- 
riaga 2000; Jaffc et al. 2001; Kinney et al. 2000; 
Bond et al. 2000; Bridle et al. 2000) on key cos- 
mological parameters. 



Furthermore, the new BOOMERANG spec- 
trum gives a value for the baryon fraction Clbh 2 
that is in excellent agreement with independent 
constraints from standard big bang nucleosynthe- 
sis (BBN), eliminating any hint of a conflict be- 
tween the BBN and the CMB-derived values for 
the baryon density (see e.g. Peebles et al. 2000). 

In this paper we support the conclusions of 
B01 by presenting two complementary analyses 
of the measured power spectrum. In Section 2 
we briefly compare these latest BOOMERANG 
results with those of BOO. In Section 3, us- 
ing "model-independent" methods similar to that 
used in Knox & Page (2000), we measure the posi- 
tions and amplitudes of the peaks in the spectrum. 
We also compute the probability distribution of 
the theoretical power spectra used in B01 given 
the measured C/s to estimate the averages and 
variances of the peaks and dips, as in Bond et al. 
(2000). 

In Section 4, we extract the distribution of cos- 
mological parameters using a different grid of the- 
oretical power spectra and a different method for 
projecting onto one and two variable likelihoods 
than that used in B01, and show good agreement 
between the two methods. The B01 C^-database 
used the cosmological parameter set {fltot, ^A, 
flbh 2 , £l c /i 2 , n s , tc, Cio}, where n s is the spec- 
tral index of primordial density fluctuations, fl c 
and f^A are the cold dark matter and vacuum en- 
ergy densities in units of the critical density, and 
Cio is an overall normalization of the power spec- 
trum. The alternate grid described here uses f2{,, 
Q m = Q c + Qt and the Hubble parameter h in- 
stead of ilfc/i 2 , Fiji 2 and Vltot — Q m + ^a, which 
are now derived parameters. 

The determination of cosmological parameters 
is affected by the presence of near-degeneracies 
among them (Efstathiou & Bond 1999). In B01 
and Lange et al. (2001), this is improved through 
the use of parameter combinations which minimize 
the effects of these degeneracies, except for the 
important one between fi tot and &a- However, 
as long as the database is sufficiently extensive 
and finely gridded, one should get nearly the same 
answer. In addition to using different parameter 
choices to generate one and two dimensional like- 
lihood functions, we use likelihood maximization 
rather than marginalization (integration) over the 
other variables. Maximization could have a dif- 
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ferent (and often more conservative) response to 
the presence of degeneracies in the model space. 
We obtain excellent agreement between the two 
treatments, and also show that when we apply our 
methods to the DASI data, we get the same results 
as Pryke et al. (2001). In Section 5 we report our 
conclusions. 

2. The Power Spectrum: Comparison 
with Previous Results 

The first results from the 1998/99 flight of 
BOOMERANG were reported one year ago by 
BOO. These included the first resolved images of 
the CMB over several percent of the sky, and a 
preliminary power spectrum based on analysis of 
a small fraction of the data. The power spectrum 
obtained by B01 incorporates a 14- fold increase in 
effective integration time, and an ~ 1.8-fold in- 
crease in sky coverage. Thus, the B01 result ob- 
tains higher precision both at low £, where the pre- 
cision is limited by sample variance, and at high 
£, where it is limited by detector noise. Moreover, 
the B01 results make use of an improved pointing 
solution that has led to qualitative improvements 
in our knowledge of the physical beam and our 
understanding of the pointing jitter, which adds 
in quadrature with the physical beam to give the 
effective angular resolution of the measurement. 

We now understand that the pointing solution 
used in BOO produced an effective beam size of 
(12.7±2.0) arcmin, rather than the (10±l)arcmin 
assumed in that analysis, due to an underestimate 
of the pointing jitter. The effect of underesti- 
mating the effective beamwidth is to suppress the 
power spectrum at small angular scales. However, 
this is not the reason preventing the detection of 
a second peak from that dataset. In fact, when 
the effects of the pointing jitter are corrected for, 
the signal-to-noise ratio of the BOO spectrum at 
t > 350 is still insufficient to detect the second 
peak, and the data are still compatible with flat 
bandpower. In Fig. 1 we plot the spectrum of BOO 
corrected for the jitter underestimate and for an 
overall gain adjustment (+20% in Ci, within the 
±20% uncertainty assigned to the absolute cali- 
bration in both BOO and B01). The corrected BOO 
spectrum is in excellent agreement with the B01 
spectrum. 

The data of BOO have been used to derive Q, to t 



in the same paper, and to derive a full set of cos- 
mological parameters in Lange et al. (2001). The 
cosmological result of BOO - that the geometry of 
the universe is nearly flat - is not altered by the 
correction, since the correction does not move sig- 
nificantly the location of the first peak. That anal- 
ysis was done assuming a " medium" prior on h and 
f2f,/i 2 . When the full parameter analysis is done, 
assuming weaker priors as in Lange et al. (2001), 
the effect of the correction is to drive fltot even 
closer to unity. With the new data set of B01, 
Q tot = lmtom usin § the e < 600 P art of thc ncw 
data, n tot = 1.02to °5 with all of it. 

The largest effect of the correction is on the 
baryon density derived in Lange et al. (2001). 
The uncorrected spectrum gives a baryon density 
tt b h 2 = 0.036±£g$j for thc "weak" prior case. The 
value derived from the B01 spectrum if only the 
data at t < 600 are used is 0.027^;^. Adding 
the data at higher I removes degeneracies and fur- 
ther reduces the value to 0.022+q;qq 3 , as described 
below. 

3. Significance and Location of the Peaks 
and Dips 

In this section we investigate the significance 
of the detection of peaks and dips in our data, 
and estimate their locations using both model- 
independent, frcqucntist methods and model- 
dependent, Bayesian methods. 

3.1. Model-Independent Analyses of Peaks 
and Dips 

3.1.1. Does a Flat Spectrum Fit the Data? 

As a first step, we answer the question, "How 
likely is it that the measured bandpowers Ct are 
just fit by a first order polynomial Cf = Ca+CbZ" ■ 
Since the first peak is evident, we limit this anal- 
ysis to the data bins centered at 450 ^ I ^ 1000. 
We compute X 2 = (C h - (%)M£(C V - where 
Mf,b' is the covariance matrix of the measured 
bandpowers Cb and is an appropriate band av- 
erage of Cg. For this exercise, we used a Gaussian 
distribution in the Cjf, but have also checked a 
lognormal in Cjf, with M suitably transformed. 
We find similar answers using these two limits of 
thc offset-lognormal distribution for the bandpow- 
ers recommended by Bond, Jaffe & Knox (2000). 
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Table 1 

Peaks and dips: location and amplitude. 









Features 


model independent 


no priors 


weak 


model independent 


no priors 


weak 


Peak 1 
Dip 1 
Peak 2 
Dip 2 
Peak 3 


213+^ 

416+12 
541+2" 

750+?° 

845+^ 


220 ± 7 
413 ±12 
539 ± 14 
688 ± 22 
825 ± 21 


221 ±5 
412 ±7 
539 ±8 
683 ± 23 
822 ±21 


040U -1100 

1850^0 
2220+18 
1550±gg 

2090+™ 


5440™ 

2480+^ 
1 530+ 470 
2170+™ 


&i4U -1030 

1590+1° 
2420+S30 

1590±Sgg 
22 70+5?o 


Dip 3 




1025 ± 24 


1024 ±26 




880+^ 


910+ 310 


Peak 4 




1139 ±24 


1138 ±24 






ii3oi«8 


Dip 4 




1328 ±31 


1324 ±33 




570+ 22 ° 


610+?« 


Peak 5 




1442 ± 30 


1439 ± 32 




680+^8 


730t£8 


Dip 5 




1661 ±37 


1660 ±36 




300lS» 


320ljg° 



Note. — Location (columns 2-4) and amplitude (columns 5-7) of peaks and dips in the power spectrum of 
the CMB measured by BOOMERANG. In column 2 we list the values measured by means of a parabolic fit 
to the bandpowers (A£ = 50). In columns 3 and 4 we list the values estimated by integrating the peak and 
dip properties of the theoretical spectra used in Lange et al. (2001) and B01 over the probability distribution 
for the database, assuming either the no prior or weak prior restrictions of those papers. In column 5,6,7 we 
do the same for the amplitude of the peaks. Note that Peaks 4 and 5, and Dips 3, 4, 5 are all outside the 
range directly measured, so they are forecasts of what is likely to emerge if the database has components 
that continue to describe the data well, as they do now. All the errors are at l-er and include the effect of 
gain and beam calibration uncertainties. 2-er confidence intervals can be significantly broader than twice the 
l-er intervals reported here, as is clear from Fig. 2. 
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Fig. 1. — The angular power spectrum of the CMB measured by BOOMERANG. The filled circles are the 
spectrum reported by B01; the up and down triangles indicate the systematic error in B01 resulting from 
the error in the beam. Three peaks are visible at multipoles ~ 210, ~ 540, ~ 840. The open squares are the 
spectrum reported by BOO, based on an analysis of ~ 7% of the data analyzed in B01. The filled squares 
present the BOO spectrum after correction for pointing jitter and scaled by 10% (la) in overall calibration 
(see text). 



We vary Ca and Cb to find the minimum of the 
X 2 , Xmin' an( i compute the probability of having a 
X 2 < Xmin- A first order polynomial is rejected in 
all cases. P( X 2 < X^-J is 95.2%, 96.5%, 94.8% for 
the multipole ranges 401-1000, 401-750, 726-1025 
respectively (using a binning of A£ = 75 for the 
bandpowers). These conclusions are robust with 
respect to variations in the location and width of 
the £-ranges, as well as for variations of the beam 
FWHM allowed by the measurement error. We 
conclude that the features measured in the spec- 
trum are statistically significant at approximately 
2a. 



3.1.2. Peak and Dip Location and Amplitude 
Likelihood Maps 

We now describe some methods we have used to 
locate the peaks and dips of the power spectrum. 
Results shown in Table 1 used a parabolic fit for Ce. 
(As described below, we have also used parabolas 
in lnC^. The use of more complex functions, see 
e.g. Knox & Page (2000), is not required by the 
data.) The criterion for robust peak/dip detec- 
tion is that the mean curvature of the parabola be 
above some threshold (e.g., some multiple of the 
rms deviation in the curvature). For the model- 
independent entries of Table 1 , for each amplitude 
C p and location of the peak £ p , we found the cur- 
vature which gave minimum x 2 (C p ,£ p ). We re- 
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Fig. 2. — A\ 2 contours for the position and ampli- 
tude of the peaks and dips in the BOOMERANG CMB 
temperature power spectrum. The crosses in the two 
upper panels correspond to the minima for the x 2 ; 
the three contours are plotted at Ax 2 = 2.3, 6.17, 11.8 
(corresponding to 68.3%, 95.4%, 99.7% confidence in- 
tervals for a gaussian likelihood in two parameters). 
The two upper panels refer to a A£ — 50 binning of 
the power spectrum. The multipole ranges used for 
the parabolic fits on peak 1, dip 1, peak 2, dip 2 and 
peak 3 are (76 - 375), (276 - 475), (401 - 750), (676 
- 875) and (726 -1025) respectively. The two lower 
panels refer to binnings with A£ = 75. While the sta- 
tistical significance of the detections is improved, the 
location of the features is less robust against the cen- 
tering of the bins: continuous and dashed lines refer 
to bin centers shifted by 35. 



peat this procedure for different data ranges. The 
contours corresponding to A\ 2 of (2.3,6.17 and 
11.8) are plotted in Fig. 2 for the ranges where 
this procedure provides a detection of both C p and 
£ p . Three peaks and two dips are clearly detected 
by this procedure, at the multipoles and with the 
amplitudes reported in the table. The curvature 
is required to be negative for peak detection, posi- 
tive for dip detection. Note that zero curvature is 
equivalent to the fiat bandpower case. Approach- 
ing this zero curvature limit is why the contours 
sometimes open up in £ p at ~ 2a and become hor- 
izontal at ~ 3cr. This is a visual reiteration of the 
point made above, that the flat Ci is rejected at 
about 2a. We have found that the location of the 
peaks is robust against variations of the gain and 
beam calibration inside the reported error inter- 
vals, 10% for the gain calibration error, producing 
a 20% error in the C p , and 1.4'(1ct) corresponding 
to 13% in the beam. The different best fit parabo- 
las for the peaks and dips of the BOOMERANG 
power spectrum are shown in Fig. 3. 

3.1.3. Fisher Matrix Approach with Sliding Bands 

We have tried a number of other model- 
independent peak/dip finding algorithms on the 
data. For the parabolic model in £ for lnC^, we 
analyzed the quadratic form in blocks of 3, 4 or 
5 bins, each bin being of width A£ — 50, de- 
termining the best fit and the variance about it 
as defined by the inverse of the likelihood curva- 
ture matrix (Fisher matrix). For this exercise, we 
adopted a lognormal distribution of the bandpow- 
ers (but checked for robustness by also assuming 
a Gaussian distribution) and marginalized over 
the beam uncertainty. We slid across the £-space 
with our fixed bin group and estimated the sig- 
nificance of the peak or dip detection by the ratio 
of the best-fit curvature to the rms deviation in 
it. For BOOMERANG, using three bins of width 
A£ = 50 gives stronger detection, but results in 
larger error bars (as estimated from the inverse 
Fisher matrix). Here, we quote our 5-bin results, 
and among these we use for each peak/dip the 
5-bin template which gives the largest ratio of 
the mean curvature to standard deviation of the 
curvature. 

We find interleaved peaks and dips at £ = 
215 ±11, 431 ±10, 522±27, 736±21 and 837±15. 
The error bars quoted correspond to the variance 
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in peak position £ p given the curvature of the 
parabola fixed at the best-fit value ( if one at- 
tempts to marginalize over the curvature, £ p loses 
localization due to contribution from fits with 
vanishing curvature, i.e. straight lines). The am- 
plitudes of the features are 

2290±!°, 1640+^° and 2210±§2o M^ 2 > corre- 
spondingly. When the 10% calibration uncertainty 
is included, the errors are similar to those in Table 
1, but of course the estimates move up and down 
together in a coherent way, so here we have chosen 
to indicate the purely statistical error bars. The 
significance of the detection as estimated by the 
distinction of the best-fit curvature from zero is 
1.7ct for the second peak and dip, and 2.2 a for 
the third peak. 

3.2. Model-Dependent Analysis of Peaks 
and Dips 

3.2.1. Ensemble- Averaging Method 

We have also estimated peak/dip positions and 
amplitudes on a large set of prescribed Ce shapes 
weighted by the probability that each shape has 
when confronted with the data. For the cases 
shown in Table 1, the prescribed shapes are the 
elements of the Ce database used in Lange et al. 
(2001) and B01. Explicitly, we compute for the 
BOOMERANG+DMR data the ensemble aver- 
ages < In £ p > and < In C p > and their variances, 

< (Aln^p) 2 >, and < (AlnC p ) 2 >, with respect 
to the product of the likelihood function for the 
parameters and their prior probability. Two cases 
are shown, for the no prior and weak prior cases 
of Lange et al. (2001). Note the good agreement 
using this method with the model-independent re- 
sults. (The results are nearly the same if we use 

< £ p > instead of < ln^ p >.) When this method 
was applied to all of the CMB data published be- 
fore 2001, including that of BOOMERANG and 
Maxima, a first peak location of 212 ± 7 was ob- 
tained (Bond et al. 2000). 

3.2.2. Forecasting Peaks and Dips 

One virtue of this procedure is that, because we 
know the spectral shape for all £, we can forecast 
where peaks and dips will lie beyond the ^-region 
we observe. Thus, Table 1 gives predictions for the 
locations and amplitudes of the subsequent (4,5) 
peaks and (3,4,5) dips. 



To test whether this works using just the data 
at hand, we restrict ourselves to using only the 
£ < 600 data and see how well the peaks be- 
yond are found. As expected, the features with 
£ < 600 are quite compatible with those in Ta- 
ble 1 derived using all of the data: for the weak 
prior case, the first peak, the first dip and the 
second peak for I < 600 are found to be at £ = 
219 ± 5, 403 ± 10 and 530 ± 23, with amplitudes 
6000t^io> 1750^0 and 2700±^° [iK 2 . The sec- 
ond dip and third peak are out of range, predicted 
to be found at I = 640 ± 49 and 789 ± 40, with 
amplitudes 2080±^ and ^^-vm f-K 2 , in good 
agreement with the £ = 683 ± 23 and 822 ± 21, 
with amplitudes 1590t|$ and 2270±|f^ fiK 2 ac- 
tually found. If we further restrict the prior prob- 
abilities, the forecasted positions can move a bit: 
e.g., we find the I < 600 cut with the weak prior 
gives iltot = 1.04 ±0.07, suggesting a constraint to 
the theoretically-motivated flat universe is reason- 
able. If we adopt as well the "large scale structure" 
prior of Lange et al. (2001), we get £ — 653 ± 45 
and 798 ±41 at amplitudes 1950^ and 29001^ 
/j,K 2 ; and the fourth peak for £ < 600 is forecasted 
to be at £ = 1107±46, similar to the £ = 1138±24 
forecast of the table. Our conclusion is that even 
if we restrict ourselves to I < 600, the forecasts are 
very good. We expect that using all of the data to 
I = 1000, with multiple peaks and dips, the fore- 
casts of the Table should be even more accurate. 

3.3. Peak and Dip Finding in the DASI 
data 

We have applied the sliding band procedure of 
Section 3.1.3 to the DASI data using the 3-band 
and 4-band sliders. For the 3-band, we find the 
first two peaks and interleaving dip at I = 202±15, 
407 ± 14 and 548 ± 10, with amplitudes 52001^, 
1590tig and 2760 ±250 ^K 2 . Again, we have not 
included the coherent 8% calibration uncertainty 
(16% in Ce) in these C p numbers. It may seems 
quite incongruous that the error bars on the peak 
and dip positions can be so small relative to the 
bin size, but we emphasize that these give the er- 
ror contours in the immediate neighborhood of the 
maximum, as described by the Fisher matrix. Just 
as in Fig. 2, the contours open up at levels lower 
than lcr, resulting in imprecise localization at the 
2(7 level. DASI's detection of the second dip is 1.6a 
in curvature by the Fisher error accounting, but it 
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Fig. 3. — Best fit parabolas for the different peaks and dips of the BOOMERANG spectrum. 



is found, at £ = 656 ± 19 and C d = l650tH° o fiK 2 . 

We have also applied the prescribed shape 
method of Section 3.2.1 to the DASI data (Halver- 
son et al. 2001). We find that the location and 
forecasts are quite compatible with those given in 
Table 1 for BOOMERANG. For example, for the 
weak prior, the first three peaks and interleaved 
two dips are at I = 216 ± 6, 401 ± 10, 523 ± 12, 
655 ±30 and 794 ±30, with amplitudes 5240±^, 

330 z, ' ±c "~'-520 

% calibration uncer- 
tainty. The forecasted third dip is at £ — 988 ± 34 
with amplitude 980^2ooV-^ 2 > compatible with the 
BOOMERANG result. 

4. Robustness of Cosmological Parame- 
ters 

The multiple peaks and dips are a strong pre- 
diction of the simplest class of adiabatic inflation- 
ary models, and more generally of models with 
passive, coherent perturbations (e.g. Albrecht et 



1570t^o, 2410^0, 1690^io an d 2480t 66i! 
These C p now include the 



al. 1996). Although the main effect giving rise 
to them is regular sound compression and rarefac- 
tion of the photon-baryon plasma at photon de- 
coupling, there are a number of influences that 
make the regularity only roughly true. Nonethe- 
less, the "catalogue" of peak and dip positions 
used to construct the table could be searched to 
find best-fitting sequences and the associated cos- 
mological parameters giving rise to them. Indeed 
we know that peaks and dips and only a few points 
in between are enough to characterize the mor- 
phology of the Cf spectra (e.g. Sigurdson & Scott 
2000). However, it is clear that it is better to work 
with the full shapes to test the theories. In this 
section, we show that the extracted cosmological 
parameters using the full shapes are robustly de- 
termined, by comparing results for the Ci database 
and Bayesian marginalization techniques used in 
Lange et al. (2001) and B01 with those obtained 
using the variables and likelihood maximization 
techniques described below and in de Bernardis et 
al. (1997), Dodelson & Knox (2000), Melchiorri et 
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al. (1999) and Balbi et al. (2000). We also test our 
methods on the DASI data set. 

4.1. Extraction of Cosmological Parame- 
ters 

The {f2 m , fife, f^A, h, n s , tc, Cio} space has pa- 
rameters sampled as follows: TO = 0.11, 1.085, 
in steps of 0.025; Q b = 0.015, 0.20, in steps 
of 0.015; n A = 0.0, ...,0.975, in steps of 0.025; 
h = 0.25, 0.95, in steps of 0.05; spectral in- 
dex of the primordial density perturbations n s = 
0.50, ...,1.50, in steps of 0.02, t c = 0.,..,0.5, in 
steps of 0.1. The overall amplitude Cio, expressed 
in units of C^f BE , is allowed to vary continu- 
ously. The theoretical models are computed us- 
ing the CMBFAST program (Seljak & Zaldarriaga 
1996), as were those used in B01. Here, we have 
ignored the role gravity waves may play. We used 
the new BOOMERANG anisotropy power spec- 
trum expressed as 18 bandpowers from I = 70 to 
t = 1050 (see B01) and we computed the likeli- 
hood for the cosmological models as exp(— x 2 /2), 
where % 2 is the quadratic form defined in section 
3.1.1. A 10% Gaussian-distributed calibration er- 
ror in the gain and a 1.4' (13%) beam uncertainty 
were included in the analysis. The COBE-DMR 
bandpowers used were those of Bond, Jaffe, and 
Knox (1998), obtained from the RADPACK dis- 
tribution (Knox 2000). For the other Ci database 
in which flk = 1 — fltot replaces d m and Q c h 2 and 
fl b h 2 replace h and the parameter grid, the 
treatment of beam and calibration uncertainty, the 
use of the offset-lognormal approximation, and the 
marginalization method, are as described in Lange 
ct al. (2001). 

As discussed in Lange et al. (2001), apart 
from the inevitable database discreteness, assum- 
ing a uniform distribution in the variables of the 
database is tantamount to adopting different rela- 
tive prior probabilities on the variables. Of course, 
this is not an issue for the maximization method, 
and even for marginalization is a very weak prior 
relative to the strong observed detections. We 
can use further "top hat" priors to mimic the 
{fl m , fib, h} database by restricting the wide cov- 
erage we have in the {flk, fi b h 2 , fl c h 2 } variables. 
We illustrate this in Fig. 4, where this mimicking 
prior is coupled to a weak cosmological prior, re- 
quiring the age of the Universe to be above 10 Gyr 
and the Hubble parameter to be in the range 0.45 



to 0.85. This is very similar to the weak prior of 
Lange et al. (2001), Jaffe et al. (2001) and B01, 
except the upper limit was 0.95: this results in 
only tiny differences in the extracted parameters. 

4.2. Cosmological Parameter Results 

Our basic results, on method testing and cos- 
mological implications, are shown in Fig. 4 and 
the associated contour plots Figs. 5 and 6. Fig. 4 
clearly shows that it does not make that much 
difference in constructing the likelihood function 
for a target cosmological variable if we marginal- 
ize (integrate) over the other cosmological param- 
eters or find the maximum likelihood value. Nor 
does it matter which database is used. Similar suc- 
cess in agreement is found with other prior choices, 
but for this paper we will restrict ourselves to the 
0.45 < h < 0.85 weak prior. When we integrate 
the distributions of Fig. 4 to get the 50% value 
and the 1-sigma errors derived from the 16% and 



84% values, we obtain Si- 



tot 



1.02 



+0.06 

-0.05 



, 1.02 



+0.05 

-0.06' 



1.04 ±0.05, for marginalization and maximization 
in the Lange et al. (2001) C^-database and the 
maximization in the alternate database, respec- 
tively; we also obtain n s = 0.961°^°,, 0.901°,;°^, 
0.90 ± 0.06, and n b h 2 = 0.022±°'°°4, 0.020^ ^, 
n mq+ 005 

U.U1»_Q QQ 4 . 

We have also applied the B01 Ci database with 
marginalization to the DASI data, and find for 
our 0.45 < h < 0.95, age > 10 Gyr, prior Q tot 



1 05+ 06 n 
i - uo -0.06' a s 



1.02l£^ and fl b h 2 



0.024 



-0.005 
-0.004' 



to be compared with the Pryke et al. (2001) values 
of 1.05±g;g|, 1.01±g;g? and 0.0221°;°°*. As men- 
tioned above, the overlap of 20% of the DASI fields 
one the sky precludes a rigorous joint analysis of 
the datasets. If you proceed anyway as an exercise 
and ignore those correlations, the results are very 
close to the BOOMERANG+DMR marginaliza- 
tion values given above, with very slightly reduced 
errors. 

We now discuss the implications of these deter- 
minations in more detail. If we knew all other vari- 
ables, the position of the first peak would allow us 
to determine the mean curvature of the universe, 
and all subsequent dips and peaks would have to 
follow in a specific set pattern. However, degenera- 
cies among cosmological parameters are present. 
In particular, there is a geometrical degeneracy 
between Q\ and fltot in the angular diameter dis- 
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Fig. 4. — Likelihood curves derived from the BOOMERANG and COBE/DMR datasets as a function of O tot , 
n s and Qbh 2 show the relative insensitivity to whether marginalization over the other variables is done (solid) 
or the maximum likelihood point in the other variables is chosen (dashed). These curves were constructed 
using a prior on the C^-database used in B01 designed to mimic the ranges in the alternate C^-database that 
uses {fib, f2 m , h} (open circles). A further weak prior, with 0.45 < h < 0.85 and age > 10 Gyr, was used. The 
high degree of correlation of n s with tq accounts for the wider distribution of the marginalized n s likelihood 
compared to that for the maximization procedure. 



tance, which allows the peak/dip pattern to be 
reproduced by different combinations of the two. 
The peak/dip heights also have near-degeneracies 
associated with them, but these can be strongly 
broken as more peaks and dips are added. Ex- 
treme examples are closed models dominated by 
baryonic dark matter, which can reproduce the 
observed position of the first peak (Griffiths et al. 
2001) but are unable to account for the observed 
I > 600 power. 

The classic Q m vs. J7a plot in Fig. 5 shows 
68.3%, 95.4% and 99.7% contours, defined to be 
where the likelihood falls to 0.32, 0.05 and 0.01 
of its peak value, as it would be for a 2D multi- 
variate Gaussian. The phase space of models in 
the fi m vs. 57a plane results in the likelihood for 
the total energy density of the universe fitot being 
skewed towards closed models (see e.g. Lange et 
al. 2001; Bond et al. 2000). This skewness is not 
important if the acceptable model space is well lo- 
calized, which can be accomplished by imposing 
prior probabilities, in particular on h. The weak 



"top hat h — 0.65 ± 0.20" prior we have adopted 
decreases the skewness to closed models (which 
is quite evident if models with very low h values 
are included, as described in B01, and helps to 
break the geometrical degeneracy, hinting at the 
presence of a cosmological constant at the level 
of ~ la. This becomes more pronounced with a 
stronger prior on h (B01). Including the recent 
supernovae data Perlmutter et al. (1997) with the 
weak prior used here, we find ft a = 0.71 ±0.11 
and Q m = 0.31^q'^2 using maximization, to be 
contrasted with the marginalization results using 
the B01 database of 0.73±[|;?£ and 0.32 ± 0.06. 

The measurement of the relative amplitude of 
the peaks in the CMB spectrum provides impor- 
tant constraints on the physical density of baryons 
fib ft 2 . In the region from I ~ 50 and up to the 
second peak, fli,h 2 is nearly degenerate with vari- 
ations in the primordial spectral index of scalar 
fluctuations n s : increasing f^ft 2 increases the ra- 
tio of the amplitudes of the first and second peaks, 
but so does decreasing n s . Beyond the second 
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Fig. 5. — Constraints in the Q m vs. fl\ plane from 
the combined BOOMERANG and COBE/DMR 
datasets, assuming the weak prior 0.45 < h < 
0.85, age > 10 Gyr. For this plot, the likelihood 
at each point is calculated by maximizing over the 
remaining 5 parameters. The shaded regions cor- 
respond to the 68.3%, 95.4% and 99.7% confidence 
levels. The CMB contours (labeled B98) are over- 
laid on the contours derived from observations of 
high redshift supernovae (labeled SN la). The line 
contours are for the combined likelihood. 

peak, however, the two effects separate: the third 
peak rises by increasing flbh 2 and lowers by de- 
creasing n s . Thus, though Fig. 6 shows that 
the two are still well-correlated, the inclusion of 
BOOMERANG data up to £ ~ 1000 has sharp- 
ened our ability to independently estimate the 
two. 

Of course, both fl^h 2 and n s are extremely im- 
portant for our understanding of the early uni- 
verse. The Qi,h 2 value can be inferred from ob- 
servations of primordial nucleides under the as- 
sumption of standard BBN scenarios. Recent ob- 
servations of primordial deuterium from quasar 
absorption line systems suggest a value tl^h 2 = 
0.020 ±0.002 at the 95% C.L. (Buries et al. 2000). 

In inflation models, the spectral index of the 



primordial fluctuations gives information about 
the shape of the primordial potential of the in- 
flaton field which drove inflation. While there is 
no fundamental constraint on this parameter, the 
simplest and least baroque models of inflation do 
give values that are just below unity. We obtain 
n s = 0.90 ± 0.09 using the maximization proce- 
dure and the n s — 0.96 ± 0.09 using the preferred 
marginalization method. 




0.00 0.01 0.02 0.03 0.04 

n b h 2 



Fig. 6. — Constraints in the fli,h 2 vs. n s 
plane from the combined BOOMERANG and 
COBE/DMR datasets. The likelihood at each 
point is computed by maximizing over the remain- 
ing five parameters. The shaded regions corre- 
spond to 68.3%, 95.4% and 99.7% confidence lev- 
els. The bound obtained from standard BBN is 
overlaid. 

Since n s and tc, the optical depth to the sur- 
face of last scattering, are also highly correlated, 
there is a noticeable difference between marginal- 
ization and maximization, as seen in Fig. 4. In- 
creasing tc suppresses the high-£ power spectrum 
by a constant factor ~ exp —2tc, but leaves the 
lowest multipoles unaffected. This is similar to 
the effect of altering the spectral index, which 
changes higher multipoles with a longer "lever 
arm" with respect to the lower £. Because the 
likelihood functions peak for nearly all models at 
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tc = 0, the maximization procedure effectively 
ignores tc > 0, whereas marginalization averages 
over the allowed values of tc, accounting for the 
increase from 0.90 to 0.96. (If we restrict ourselves 
to tc = 0, we get 0.89^Z q| with maximization, 
0.92+;^ with marginalization. ) 

The values of the spectral index, the curvature 
and the cosmological constant affect the shape 
and the amplitude of the power spectrum of the 
matter distribution. The bandpower erf, giving 
the variance in (linear) matter fluctuations aver- 
aged in 8h~ 1 Mpc spheres, is often used to char- 
acterize the results of many large scale struc- 
ture observations. In Lange et al. (2001), Jaffe 
et al. (2001) and B01, CT 8 f^ 56 =0.55+;o2;±'.S8 was 
adopted, where the combination of Gaussian and 
top hat error bars were used to generate a wide 
distribution "weak" large scale structure prior to 
be imposed upon the CMB data. This compares 
with er 8 ^ 5 ~ 0.50 ±0.05 estimated by Pen (1998) 
and a best fit value ergSl^ 47 ~ 0.56 from Viana 
& Liddle (1999). We have computed cr 8 for each 
model in our database using the matter fluctua- 
tions power spectrum. From the likelihood analy- 
sis in the space (as, TO ), we find the 1-sigma range 
0.5 < cr 8 < 0.8 for 0.3 < Q m < 0.7. For the de- 
tailed application of the weak large scale structure 
prior of Lange et al. (2001), which also included 
a constraint on the shape of the power spectrum, 
see B01. 

5. Discussion and Conclusions 

We have analyzed the most recent results from 
BOOMERANG, and shown that there are a series 
of 5 features, 3 peaks with 2 interleaved dips, each 
of which is detected using a cosmological model in- 
dependent method at ~ 95% or better confidence. 
Because BOOMERANG is able to achieve resolu- 
tion in ^-space of A£ ~ 50, with only ~ 10% (auto- 
correlation between bins, the data determine the 
positions as well as the amplitudes of each of the 
features with reasonable precision. 

The positions and amplitudes of the features 
are consistent with the results of the DASI exper- 
iment. A direct, cosmological model- independent 
comparison of the two experiments is made diffi- 
cult by the lower £-space resolution of the DASI 
experiment. However, the two results can be ac- 
curately compared via a model-dependent analysis 



that takes millions of Ce shapes with very differ- 
ent peak/dip amplitudes and locations, weighted 
their peak and dip locations with the probability of 
the (normalized) shapes, and finds well- localized 
£ p and C p for each of the peaks and dips. The 
agreement between BOOMERANG and DASI us- 
ing this method is very good, and completely con- 
sistent with the model-independent results. It will 
of course be of great interest whether the forecasts 
allowed by this model-dependent technique for the 
positions and amplitudes of the next peaks and 
dips will be borne out by even higher resolution 
observations. 

The natural interpretation of this much-antici- 
pated sequence of peaks and dips in d is that 
we are seeing phase-coherent pressure waves in 
the photon-baryon fluid at photon decoupling 
at redshift <~ 1100, expected in adiabatic mod- 
els of structure formation, and which therefore 
the simplest inflation models give. This is cer- 
tainly the most economical interpretation, espe- 
cially since the locations within this paradigm 
correspond to widely-anticipated cosmological pa- 
rameter choices, namely VL to t ~ 1, ~ 1/3, 
Q, B h 2 w 0.02, with a nearly scale invariant spec- 
trum, n s <~ 1. We have shown how well deter- 
mined and robust these parameter values are to 
changes in the Ci database, in using likelihood 
maximization or Bayesian marginalization, and to 
using either the BOOMERANG or DASI C £ s in 
conjunction with the DMR Ci&. 

Of course this does not clinch the case: the de- 
rived cosmological parameters could be quite dif- 
ferent if we allowed ourselves much further free- 
dom beyond a single slope to characterize the pri- 
mordial spectrum. It could even be that the peaks 
and dips reflect early universe structures rather 
than sound wave structures. It is also possible 
that isocurvature modes with artful enough initial 
condition choices could mimic the simple adiabatic 
case (e.g. Turok 1996). The detection of the asso- 
ciated polarization peaks and dips would rule out 
these more exotic possibilities. Although there are 
no glaring anomalies between theory and data at 
this stage requiring a revisit of basic CMB assump- 
tions, the improved precision in this £-range that 
more sky coverage will give, and the extension to 
higher t that other CMB experiments will give, 
could well reveal that our forecasts are wrong. 
Even with the current data, there is certainly room 
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for other cosmological parameters not treated in 
our minimalist inflation databases, and such ex- 
plorations are underway. 
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